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ABSTRACT 


An  impco/ei  an  I  simplified  model  foe  the 
Heave-Pitch  dynamics  of  the  XR-3  CA3  Surface  Effect 
Ship  is  davaLoped  for  constant  speai  operation.  The 
nonlinear  equations  of  notion  ara  linaarizei  aoout  the 
steady-state  operating  point.  rime-iomain  validation 
is  aocomplishai  by  comparison  with  taa  6  DOF  nodal,  A 
signal  flow  graph  of  the  craft  dynamics  is  ievaloped 
and  Mason's  3ain  Rule  used  co  determine  the 
characteristic  S-Oomain  polynomials  describing  the 
craft's  vertical  pLana  dynamic  baha'/ior,  »ith  Bode 
Plot  analysis  included.  Conclusions  are  irawn  and 
racommandations  for  further  study  are  made. 
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INTRODUCTION 


In  1973  Leo  and  Boncal  [1]  converted  the  Oceanics  Inc.  6 
DOF  Loads  and  Motions  digital  computer  simulation  program 
for  the  100-B  Surface  Effect  Ship  to  reflect  the  physical 
characteristics  of  the  XR-3  craft.  Since  that  time  there 
has  been  a  continuous  effort  at  the  NPS  to  update  and 
improve  on  this  simulation  [2,3,4,5]  . 

One  of  the  major  Irawbacks  of  that  program  which 
precludes  regular  ass  of  it  for  sea  state  operations 
simulation  is  the  excessive  computational  and  turn-around 
times  required  [6]  .  ks  it  is  presently  iesired  to  use  the  6 
DOF  simulation  for  design  and  evaluation  studies  of  the 
proposed  3K  Ton  5SS  control  System  in  varying  sea  state 
conditions,  it  is  important  to  significantly  reduce  the 
simulation's  computational  time. 

One  area  where  the  possibility  exists  to  achieve  a  major 
reduction  in  computational  time  is  tie  manner  in  which 
sidewall  forces  and  moaents  are  calculated.  The  present 
form  of  the  6  DOF  progran  uses  a  large  number  of  discrete 
sections  (28  per  sidewall)  with  an  off-line  table  look-up 
scheme  to  determine  the  forces  generated  in  each  discrete 
section.  The  simple  model  derived  by  Molntyre  Zll  uses  two 
sections  per  sidewall  witi  very  simplified  geometry  that 
permits  on-Line  computation  of  the  forces  and  aonents.  The 
result  of  eventual  incorporation  of  this  sidewalL  geometry 
and  forces  calculation  into  the  6  DOF  program  should  be  a 
significant,  reduction  of  computation  and  turn-around  time. 

Another   area   of   great   interest   in  SES  Siuulation  is 
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gaining  insight  into  the  dynamic  behavior  of  ttie  craft's 
vertical  motion  of  pitch  and  heave  due  to  pressare  induced 
lift  forces  generated  in  the  plenum  and  seals.  As  stated  in 
ref  C8]  as  recently  as  August,  1978, 

••The  overall  precision  of  the  analysis  is  satisfactory 
for  predicting  stability  and  maneuverability.  However, 
several  areas  present  ?ery  complex  or  difficult  physical 
phenomena  to  analyze  and  further  effort  should  be 
invested  to  improve  understanding  of  these.  '  ' 

The  first  area  for  continuing  investigation  reconnead ed  was; 

•'A  oetter  determination  of  pitch-heave  dynamic 
characteristics.  '  • 

An  improved  knowledge  of  the  dynamic  characteristics 
brougat  about  by  the  pressure  induced  forces  and  moments 
could  lead  to  a  reduction  of  the  computation  time 
required  for  craft  simulation.  Ihese  pressure  effects 
introduce  a  stiffening  in  the  integration  of  sea  state 
induced  oscillations  and  a  better  understaiding  of 
processes  could  lead  to  the  removal  of  this  stiffness 
and  thus  reduced  CP0  time. 

In  this  work,  both  of  these  areas  are  considered.  First 
of  all  the  various  craft  physical  characteristics  are 
modeled,  resulting  in  changes  and  improvements  to  the 
initial  efforts  by  Mclntyre  C7J  .  These  equations  are 
linearized  about  a  steady-state  operating  point  using 
the  Taylor  Series  expansion.  The  linearized  equations 
are  cast  into  a  State  Space  format  and  simulated  on  the 
digital  computer.  Comparison  is  made  between  the 
time-domain  results  of  the  linear  model  and  tie  6  DOF 
model.  The  State  Space  representation  is  used  to 
generate  a  signal  flow   graph.    Mason's   Gain   Rule   is 


11 


applied  and  the  S-domain  roots  of  the  craft's 
characteristic  equations  are  located  and  it  is  shown 
that  soma  simplification  and  redaction  of  system  order 
is  possible.  The  resultant  S-domain  polynomials  are 
used  to  generate  Bode  Plots  to  graphically  display  the 
craft's  simplified  frequency  response. 

Conclusions  are  dra*n  and  recommendations  are  made 
concerning  the  validty  of  the  model  and  areas  for  future 
study. 

Several  procedural  notes  are  in  order.  Because  all 
equations  are  to  be  solved  digitally  using  the  FORTRAN 
computer  language  and  because  digital  computer  text 
processing  facilities  were  used  to  generate  this 
document,  all  formulas  will  be  presented  in  FORTRAN 
format.  In  the  thesis  body  reference  to  the  ''nonlinear 
model'1,  ''six  degcee-of  -freedom  model' *  or  ''6  DOF'' 
model  all  mean  the  Loads  and  Motions  program  as  adapted 
to  the  XR-3  test  craft  it  the  Naval  Postgraduate  School 
[3]  .  Reference  to  the  linear  model  refers  to  the 
linearized  model  of  the  same  craft  as  developed   herein. 
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II.   NONLINEAR  XR-3  MODEL  HEAVE  AND  PITCH  EQUATIONS 


A.   BACKGROUND 


The  original  development  work  on  the  simplified 
two-degrea-of  -freedom  model  was  dona  by  Hclntyre  CB  .. 
Several  weaknesses  in  the  original  model  wera  noted  by 
Mclntyre  in  that  work. 

The  major  portion  of  the  work  presented  here  was 
involved  in  close  examination  of  all  aspects  of  the 
original  model,  comparison  of  the  modeling  techniques 
used  therin  with  other  methods  and  remodeling  of  several 
important  craft  physical,  characteristics. 

One  major  area  of  investigation  was  the  lack  of  pitch 

damping   in   the   original   model,  which  resulted  in  the 

addition  of  an  added-nass  effect   to  '  the   pitca  moment 
equations . 

Additional  changes  involve  a  revision  of  the  force 
aquations  for  the  rear  seal,  major  changes  in  the  way 
the  plenum  chancer  is  modeled  (also  added  to  the  loads 
and  motion  6-dof  modal  C2H  )  .  and  a  more  rigorous 
development  of  the  planing  force  equation  of  [7]  . 

There  are  saveral  areas  of  tneory  concerning  CAB  5ES 
modeling  where  empirical  data  is  lacking  due  to 
difficulties  with  physical  measurements.  One  such  area 
inviting  closer  inspection  was  dynamic  variation  of   the 
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effective  center  of  pressure  foe  the  plenum  gauge 
pressure  lift  effect.  I  possible  approach  is  suggested 
herein  and  included  in  the  model. 


B.   ASSUMPTIONS 


The  following  assumptions  are  made  for  this  constant 
speed,  pitch-heave  model: 

1.  Planing  forces  are  concentrated  it  the  centroid  of 
the  keel  longitudinal  cross  section. 

2.  Total  seal  force  is  lumped  into  a  pressure 
differential  acting  on  the  seal  face  area  in  contact 
with  the  water  passing  beneath  the  seal. 

3.  Aerodynamic  faeces  on  the  forward  seal  are 
disregarded. 

4.  Dynamic  variations  in  moment  lever  arms  are 
disregarded  for  all  moaeats  except  seals  and  the  moment 
due  to  plenum  guage  pressure. 

5.  Dynamic  variations  in  effective  plenum  roof  area 
will  be  included. 

6.  Any  lift  effects  lue  to  aerodynamic  forces  are 
lumped  in  with  the  planing  forces  (for  reasons  to  be 
developed  in  Chapt  IV.)  The  moments  are  incLuded  in 
plenum  gauge  pressure  moments  . 

7.  Constant  seal  leakage  area  is  assumed. 

3.   Since   this   aodeL  is  for  constant  surge  velocity, 
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drag   and   thrast   forces   are   considered   to    be    in 
equilibrium  and  in  balance  for  moments  about  the  1    axis. 


C.   COORDINATE  SYSTEM 


(see  fig.  1  ) 

1.  The  origin  of  the  cartesian  coordinate  system  is 
located  at  the  call-water  line,  longitudinally  and 
transversly  coincident  with  the  craft  center  of  gravity. 

2.  The  X  axis  is  the  longitudinal  coordinate  defined 
as  positive  as  X  increases  towards  the  bow  of  tie  craft. 

2.  The  Y  axis  is  the  transverse  coordinate  increasing 
positively  towards  the  starboard  side. 

4.  The  Z  axis  is  the  vertical  axis,  increasing 
positively  in  a  downwari  direction. 

5.  ZS  is  the  vertical.  distance  from  the  center  of 
gravity  to  the  keel. 

5.  Z  is  ~he  vertical  distance  from  the  calm  water  line 
to  the  vertical  center  of    gravity. 

7.  THETA  is  the  pitch  angle  defined  as  positive  as  the 
bow  pitches  up. 

8.  Positive  pitch  moments  are  defined  as  moments  which 
tend  ;o  cause  a  positive  increase  in  the  angle  THETA. 


D.   SIDEWAH.  MODELING 
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As  it  was  desired  to  keep  all  physical  dimensions 
consistant  with  the  six  degree  of  freedom  modal,  scale 
model  drawings  ware  developed  with  the  data  from  the 
sidewall  subroutine  of  that  modal.  All  sidewall 
dimensions  usad  in  this  model  wera  taken  from  those 
drawings. 

The  sidewalls  are  modaled  in  two  sections,  ona  forward 
of  the  longitudinal  Z.Z.  and  the  other  aft.  Each 
sidewall  is  considered  to  be  of  constant  width  at  the 
water  line  and  the  ksel,  with  width  increasing  in  the 
vertical  direction  according  to  an  average  deadrise 
angla. 

The  average  forward  sidewall  area  is  found  to  be: 

(see  fig.   2   and  3b  ) 

AV2RAGS  WIDTH  =  LDB&H1  /  (2*TAN  (DR  1)  )  +  WS1  0 

(II-1) 

whare  LDBAR1  is  draft  at  forward  centroid,  DR1  = 
average  deadrise  angle  forward  and  MS  10  is  the  average 
keel  width  of  the  forward  section.  Therefore  submerged 
volume  of  a  forward  sidewall  is  given  by: 

VOL  =  L1*  (LDBAR1/(2*TAN (DR1) )  +WS10)  *LDBAR1 

(II-2) 

where  L1  is  the  length  of  a  forward  buoyant  section. 
Following  the  same  procedures  the  aftar  sidewall  width 
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is: 


(see  fig.   2   and   3c  ) 


AVERAGE  WIDTH  =  LDBAR2/ (2*TAN (DR2)  )  + WS20 


(II-3) 


where  DR2  =  average  deadrise  angle  aft 
LD3AR2  =  draft  at  after  buoyant  centroii 
HS20  =  average  keel  width  aft 

Submerged   volume   of   an   after   sidewall    is    then 
determined  to  be; 


VOL  =  L2*  (WS20  +  LD3AR2/(2*TAN(DR2)  )  )  *LD3AR2 

(II-1*) 

where  L2  is  the  length  of.    the  after  sidewall. 


E.   SEAL  MODELING 


A  major  change  in  Hclntyre's  model  was  made  by 
remodeling  the  stern  saal.  In  the  original  model  a  large 
negative  moment  was  required  to  place  the  pitch 
equations  in  steady  state  equibibrium.  Investigation  of 
seal  forces  as  presented  in  [5Q  revealed  a  large 
discrepancy  between  the  after  seal  forces  generated  in 
Ilclntyre's  model  aad  calculations  based  on  the 
information  in  EQ  .  This  discrepancy  was  further 
supported  by  actual  seal  load  lata  presented  by  Layton 
in  DCfl  •  All  seal  dimensions  used  in  this  model  were 
taken  from  data  presented  in  D  "Q  .  Seal  hinge  locations 
were   taken    from    the    data    used    in    G]  . 
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Seal  lift  force  forward  is  generated  by  planum  gauge 
pressure  acting  downward  on  the  seal  face  area  in 
contact  with  the  water  passing  underneath  the  seal. 
This  downward  force  results  in  a  lift  force  being 
generated  as  the  downward  movement  of  the  seal  displaces 
the  water  beneath  it. 

The  seal  is  considered  to  have  a  straight  forward  face 
which  runs  from  the  hinge  point  at  the  plenum  roof  to 
the  oottom  of  the  kssl  at  an  angle  determined  to  be 
best-fit  from  the  data  in  ref.  D  "0  .  The  wetted  surface 
of  the  seal  is  considered  to  la/  parallel  to  and  on  the 
water  at  the  point  where  the  seal  face  intercepts  the 
calm  water  line,  as  a  result  forward  seal  lift  equations 
are : 

ASEAL1  =  »IDTH*XSEAL1 

(II-5) 


XSSAL1  =  (LD-L3*TAM  (TH3TA)  )/SIN(31) 

(II-6) 

LD  =  draft  at  the  C.G 

L3  =  distance  from  C.G  to  point  of  contact  for  forward 

seal  face. 

LD   -    L3*TAN(TKETA)     =   draft    at    forward   seal    face    contact 

point 

31    degrees   is    the    interior    angle    of    the    seal. 

The      after   seal    wetted    area    is   determined   in    a    similar 
manner    to    be    : 

A5EAL2    =    WIDTH*XSEAL2 
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(TI-7) 

XSEAL2    =     (LD+L4*TAN  (THETA)  ) /SIN(32) 

(II-3) 

In  this  modal  forward  seal  pressure  is  considered 
equal  to  plenum  pressure.  As  a  result  forward  seal  lift 
is    : 

HSEALF    =    P3BAR*ASEAL1 

(H-9) 

Since  plenum  guage  pressure  tends  to  force  the  after 
seal  upward,  resulting  in  venting  of  the  planum  air 
mass,  the  rear  seal  is  operated  at  a  pressure  slightly 
above  that  of  the  plenu.11  to  reduce  tha  venting  tendency 
of  the  after  seal.  As  a  result  the  after  seal  requires  a 
modification  in  moialing.  In  flclntyre's  work  the 
pressure  differential  (PSE  AL-PBBAR)  only  was  considered 
in  the  lift  force.  However,  as  previously  stated,  there 
were  large  discrepancies  in  computed  loads  as  compared 
to  measured  data  from  Cl  OQ  .  Therefore  in  this  model  , 
the  seal  working  pressure  is  consiiared  to  ba  ?33AR  plus 
differential  pressure.  As  a  result  after  seal  Lift  is 
given  by: 

HSEALA  =  (PBBAR  +  PDIFF)  *A5EAL2 

(11-10) 
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PLENUM  MODELING 


In  the  work  presented  by  Mclatyre  the  plenum  was 
considered  a  box  with  the  only  dynamic  dimensional 
variation  being  in  the  Z  direction.  However  (as  can  be 
seen  in  fig.  1  ) ,  although  the  forward  edge  can  be 
assumed  to  be  reasonably  constant  in  the  Y-Z  plane,  the 
forward  face  of  the  rear  seal  slopes  rather  steeply. 
Therefore  the  plenum  length  (in  the  X  direction)  varies 
rapidly  as  draft  varies.  As  a  result  the  affective 
plenum  roof  area  consists  of  that  area  from  the  after 
edge  of  the  forward  seal  to  the  forward  hinge  point  of 
the  rear  seal,  plus  the  vertical  projection  of  the 
unwetted  forward  face  of  the  raar  seal.  As  an 
additional  consequence  it  can  be  seen  that  the  effective 
plenum  center  of  pressure  varies  fore  and  aft  as  draft 
varies  .  In  fact  as  draft  increases  the  -.P.  moves 
forward. 

The   following  variable  names  will  be  used  to  describe 
the  plenum  chamber  (see  fig.  4.) 

LPS   =   length   of  plenun  roof  from  the  rear  eige  of  the 

forward  seal  to  the  hinge  point  of  the  after  seal. 

LPWL   =  length  of  plenum  at  the  keel  from  the  after  edge 

of  the  forward  seals  to  the  transom. 

LP  =  length  of  plenum  at  the  waterline. 

3U3HGT  =  vertical  distance  from  the  plenum  roof  to  the 

bottom  of  the  keel  at  the  center  of  gravity. 

A3  =  plenum  effective  roof  area. 

Using  the  previous  dafinitions,  the  plenum  roof  length 
and  area  is: 
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LP    =    LP  ML-  (LPWL-LPR)  *  (LD/BUBHGT) 

(II-11a) 

and   effective   plenum   roof    area    is 


AB    =    MIDTH*LP 

(II- Tib) 

with  the  result  that  lift  due  to  plenum   gauge   pressure 
is: 

HPRES  =  -A3*PBBAR 

(II-11C) 

As  a  comment  on  the  complexity  of  the  lift  forces  of 
this  model  note  that  for  constant  plenum  pressure 
increased  draft  results  in  decreased  plenum  lift. 

G.   CRAFT  DYNAMICS 


1 •   Plenum  Air  Mass  ia i  Pressure  Equations 

For  the  purposes  of  this  model  the  rsar  seal 
leakage  area  is  considered  constant  since  to  assume 
otherwise  would  lead  to  modeling  a  very  difficult  and 
poorly  understood  dynamic  process.  In  addition  the 
plenum  pressure  VS  air  nass  process  is  assumed  to  be 
adiabatic  in  nature.  Fan  characteristics  are  those  used 
in  the  constant  air  leakage  heave  only  model  of  Gerba 
and  Thaler  in  D  <G  .  kir  flow  into  the  plenum  is  taken 
as  positive. 
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Therefore  the  airflow  into  the  chamber  is  given 
by: 

Qin  =  N*(QlO-Kq*PBBAR) 

(11-12) 

Where  N  is  the  nuaber  of  fans  supplying  air  to  the 
plenum,  QIO  is  the  equilibrium  air  flow  rate  for  the  fan 
with  zero  gauge  pressure  and  Kg  is  the  slope  of  the 
fans'  airflow  VS  gauge  pressure  curve  . 

Qout  is  airflow  out  of  the  plenum  and  is 
determined  oy  the  seal  leakage  equation: 

Qout  =  Cn*ai*(2*PBBAB/HH0A)**1/2 

(11-13) 

Where  Ca  =  seal  leakage  coefficient,  Al  =  seal  leakage 
area,  and  RHOA  =  mass  density  of  air.  Therefore  mass 
flow   rate    is: 

MflDOT    =    RHOA* (Qin-Qout) 

(11-14) 

The  plenum  air  dynamics  are  considered  in 
equilibrium  when  Qin  =  2out 

Absolute  Plenum  Pressure  is  determined  from  the 
adiabatic  gas  law  to  be: 

Pb  =  Pa*(MB/  (RH0A*PB3AR)  )  **GAMMA 
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(11-15) 

and, 

PB3AR    =   Pb-Pa 

(II-15a) 

Where  Pa  =  ambient  atmospheric  pressure  <LBF/FT**2  GAtfMA 
=  adiabatic  process  coefficient,  Pb  =  plenum  absolute 
pressure,    and    PBBAR    =   planum    gauge    pressure. 

2 •       Buoyant   Forces 

Buoyant  force  is  calculated  by  the  submerged 
volume  multiplied  by  the  mass  density  of  the  fluid  times 
the   gravitational   constant. 

Therefore    forward    buoyant    force   is: 

HBF    =    -2*RH0*G*V0L 

(11-16) 

Where  volume  is  given  by  aqn.  II-2,  RHQ  is  the 
mass  density  of  water  (slugs/f  t**3 )  and  3  is  the 
gravitational  constant.  The  factor  of  2  is  introduced 
to  account  for  2  sidewalls.  The  total  expression  for 
HBF  is: 

HBF  =  -2*RH0*G*L1*LDBAR1*  (LDBAR1 / (2*TAM  (DR1) )  +WS01) 

(11-17) 


Following  the  same  pattern, 
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HBA  =  -2*RH0*G*L2*LDBAR2*  (LDBAR2/ (2*TAN  (DR2) )  =W510) 

(11-18) 

3.   Seal  Forces 


be: 


Forward  seal  force  is  determined  from  eqn.  II-9  to 


HSF  =  -P3BAR*WIDTH* (LD-L3*TAN  (THETA)  ) /SIN  (31) 


(11-19) 


where    HSF    is   the   heave   force   due    to   forward   seal. 
In   a    similar   manner    HSA    is: 

HSA    =     (PBBAR  +  PDIFF)  *WTDIH*  (LD  +  L4*TAN  (THETA)/SIN  (32) 

(11-20) 


*•      Planing  H2££§£ 

Mclntyra  introduced  the  planing  force  phenomena  in 
his  model  to  account  for  the  proven  fact  that  the  CAB 
SES  tends  to  decrease  it's  draft  and  pitch  down  as  speed 
increases.  His  development  of  the  planing  equation  was 
not  well  explained  aid  the  basis  for  selection  of  the 
planing  coefficient  was  not  clearly  stated.  Raidel  D 2] 
did  a  special  investigation  of  the  planing  force 
phenomena.  His  results  were  based  on  flat  plate  planing 
effects  to  approximate  the  flat  bottom  of  the  craft 
sidewalls. 
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Since  the  water  is  in  contact  with  the  vertical 
portions  of  the  sidewalls  it  was  felt  a  more  accurate 
representation  was  neeced.  Barnaby  D3]  states  that 
water  contact  on  the  vertical  sides  of  a  planing  hull 
rapidly  reduces  planing  lift  advantage  since  one  of  the 
advantages  of  a  planing  hull  is  reduced  drag  due  to 
wetted  surface.  planing  lift  results  in  a  speed 
sensitive  tradeoff  between  static  lift  (buoyant  affects) 
and  dynamic  lift  (planing  effect)  that  is  directly 
proportional  to  the  angla  of  attack  (for  small  planing 
angles)  and  the  sguare  of  the  velocity.  Therefore  the 
planing  force  is  of  the  form: 

HPLAN  =  K*V**2*THETA 

From   Reidel's    work, 

HPLAN    =    0. 5*AREA*V**2*PI*SIN(THErA) 

(11-2  1) 
for    small   angla    approximations   this    simplifies    to: 

HPLAN    =    RHO*AREA*V**2*PI*THETA 

Which  is  the  planing  focce  for  a  flat  plate.  Barnaby 
gives  no  simple  formula  for  the  losses  due  to  Immersion 
so  we  are  left  with  an  arbitrary  loss  coefficient, 
PLCOEF  such  that: 

HPLAN  =  (PLCOEF) *2*RHO* AR5A*V**2*PI*THST A 

(11-22) 

The  added  factor  of  2  accounts  for  2  sidewalls. 

The   planing  area  is  taken  from  the  scale  drawings 
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developed  from  sidewall  data,  and  the  force  is  assumed 
to  act  at  the  centroii  of  the  flat  keel  surface,  which 
is  consistant  with  ReideL's  work.  The  loss  coafficient 
will  be  developed  in  the  section  on  aohiaving 
equilibrium. 

5«       Plesum   Gauqe   Press  are   Lift 

Combining    equations    II-11a    ,II-11b      ,      and      II-11c 
results   in: 

HPRES    =    -PBBAR*WIDTH*(L?WL- (LPWL-LPR)  * (LD/BUBHST)  ) 

(11-23) 

6.   Total  Lift  Forces 


Lift  Total  =  HPRES  +  HBP  +  HBA  +  HSP  +  HSA  *■    HPLAN 

which   must   exactly  equal  craft  weight  (W)  for  the  hull 
to  be  in  vertical  eguilibrium,  or; 


W  ••-  HPRES  +  HBF  +  HBA  +  HSF  +  HSA  +  HPLAN  =  0 

(11-24) 
defines  equilibrium  condition  in  heave  motion. 

H.   CRAFT  PITCH  MOMENT  EQUATIONS 


t •   Buoyant  Moments 
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Forward  buoyant  moment  is  the  forward  buoyant  lift 
times  the  effective  lever  arm,  or; 

PBF  =  -HBF*  (Lever  arm  forward) 

PBF  =  -HBF*L5 

(11-25) 

Note  that  the  negative  sign  results  from  the 
definition  of  positive  moments  and  the  fact  that  upward 
forces  are  negative  by  lefinition. 

PBA  =  HBA*LEVER  ARM 

P3A  =  SBA*L6 

(11-26) 
using  the  same  reasoning,  after  buoyant  moment  is: 

PB  total  =  PBF  *    PBA 

(H-27) 

2  •   S  ea  1  Mojn  en  ts 

Since  the  seals  ace  relatively  distant  from  the 
C.G.  dynamic  variations  in  the  lever  arms  are  considered 
in  seal  moment  equations  and  are  developed  nere. 
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The  seal  forca  is  modalad  to  act  at  the 
longitudinal  center  oft ha  seal  wattad  length.  Therefore 
the  effective  laver  an  for  the  seal  moment  is: 

PLSF  =  L3  -  XSEAL1/2 

whera  PLSF  =  moment  lever  arm  for  seal   forward.    As  a 
result : 

PSF  =  -HSF*PLSF 


PSF  =  -HSF*(L3-(XSEAL1/2)  ) 


(11-28) 


In   a    similar    lanaar 
PLSA    =    L4    +    XSEAL2/2 

PSA    =    tiSA*?L5A 

PSA    =    HSA*  (L4  +  XSEAL2/2) 

(11-29) 

PS  TOTAL  =  PSF+PSA 

Note  that  the  dynamic  pitch  sensitivity  to  moment 
arm  length  is  included  in  the  XSEAL  tarns  as  dafined  in 
aquations  II-5  and  II-3. 
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3 •   Planing  Moment 

For  the  planing  force  developed  in  equation  11-22 
the  planing  force  is  concentrated  at  the  centroid  of  the 
flat  keel  which  results  in: 


PPLAN  =  L7*HPLAN 


(11-30) 


*«   Plenum  Pressure  ioment 

First  the  dynamic  variation  in  effective  center  of 
pressure  will  be  developed,  (see  fig.  (4)) 

plenum  length  is  taken  from  equation  II-11a  to  be 

LP  =  LPWL-  (LPWL-LP3)  *  (LD /B UBHGT) 
Defining  LDIFF  to  be 


LDIFF  =  LPWL-LPS 


and, 


(H-31) 


XCP    =    XCP0+  (LDIFF*LD/3UBHGT)  /2 


where 


(11-32) 
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XCPO    =    LPWL/2-XCG 

(H-33) 

and      XCG      is      the      distance   from   the    transom   to    the    C.G. 
(see   fig.       4    .)       Finally, 

PPRES    =   -HPRES*XCP 

(11-34) 


5.   £it£h  Daaginq  foment 

Any  frictional  damping  terms  dependent  upon  the 
square  of  pitch  rate  (i.e.  tangential  velocity  about  the 
I  axis)  disappear  when  the  Taylor  Series  expansion  is 
carried  about  the  zero  rate  ejuilibrium  point.  Since 
for  small  signal  approximation  the  model  should  reflect 
actual  craft  characteristics,  it  was  hoped  that  the 
changes  in  seal  and  plenum  modeling  would  reveal  the 
missing  damping  characteristics.  ks  it  turnei  out  this 
was  not  the  case.  Examination  of  the  hydcodynamic 
equations  in  [93  revealed  an  added  mass  effect  dependant 
upon  a  cross  velocity  term  between  linear  X  direction 
velocity  U,  pitch  rate  Q,  and  addad-mass  at  the  stern 
A33S  due  to  slender  body  theory  applied  to  the  sidewalls 
.  This  moment  is: 

FP  =  -A33S*XLSS**2*U*TaETADOT 

(11-35) 

Where  XL5S  is  the  distance  from  the  C.G.   to  the   stern, 
and  THSTADOT  is  pitch  cate. 
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A33S  =  YSS**2*RH0*PI/8 


(11-36) 


whara  YSS  is  the  cross  section  width  of  the  side 
wall  at  tha  stern  wateriine.  Therefore, 

PPDAMP  =  -2*XLSS**2*A33S*U*THETAD0T 

=  -2*PI*XL5S**2*ILSS**2*RHO*0*THSID0T  /8 

(11-37) 

The  factor   of   2   is   introduced   to   account   for   two 
siiewalls. 

6-   £2ta.i  Bqaents 

• 

Since  tha  moments  are  all  rafarancad  to  the  craft 
center  of  gravity,  the  craft  weijht  contributes  zaro 
Moment  in  that  frame  of  refaranca.  Therefora  total 
moments  are: 

PSF  +  PSA  +  PBF  *  PBA  *■    PRE5  +  PPLAN  +  PDAMP  =  0 

(11-38) 


I.   HEAVE  ACCELERATION  EQUATIONS. 


From   Newton's  law,  acceleration  and  force  ara  related 
by: 
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FORCE  =  MASS*ACCELSRAriON 

For   the   linearized   model   we   ara   interestal   in  the 
acceleration  on  the  craft  and  as  a  result: 


ACCELERATION  =  FORCE/  MASS  and  therefore  the  heave 
equation  of  motion  becomes: 

ZDDOT  =  (HBF+HBA+HSA+HPLAN+HPRES) /MASS 

(11-39) 

where   ZDDOT   is  the  second  derivative  of  Z  witi  respect 
to  time. 

J.   PITCH  ACCELERATION  EQUATIONS 

Angular  acceleration  is: 

THETADDOT  =  PTOTAL/IYY 

Where  IYY  is  tha  craft  inartial  momant  about  tae  I    axis. 
As  a  result  total  anguLar  acceleration  is; 

THETADDOT  =  (PBF *P3A +PSF  +  PS A+PPLAN  +  PPR ES+PDAMP)  /IYY 

(11-40) 

This  completes  the  derivation  of  all  nonlinear 
aquations  naadad  to  describe  the  two-degree-o£- freedom 
CAB  SES.  In  the  following  chapter  thesa  lonlinear 
equations  will  be  linearized  using  the  Taylor  Series 
expansion  about  a  steady  state  operating  poiat  which 
will   result   in   the  equations   naeded   to   generate  a 
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dynamic   digital-computec    simulation    of    craft      notion      in 
the   time-domain. 
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III.   LINEARIZATION  OF  THE  HEAVE- PITCH  12U|riONS^ 


A.   TAYLOR  SERIES  EXPANSION 

- 

For  small  perturbations  about  a  steady  state  operating 
point  the  nonlinear  function  Z  can  be  approximated  by: 

Z  =  Z  (X,Y) 

Z(0)+dZ  =  Z(X(0)  ,1(0)  )  f  (dZ/dX)  *dZ*  (dZ/dY)  *dl 

where  OZ/dX)  and  (dZ/dY)  are  the  partial  derivatives 
of  the  function  Z  evaluated  at  the  steady  state 
operating  point  X(0)  ,Y(0)  .  Cancelling  the  steady  state 
values  from  both  sides  of  the  equation  yields  the 
differential  model  equations: 

dZ  =  (&Z/dX)  *dX+  (dZ/dY)  *dY 
and 

Z  =  Z  (0)  +dZ 

where  dZ,dX,  and  dY  ace  the  differential  variables  for 
Z,  X,  Y  respectively,  QZ/<3X)  and  (4Z/3Y)  are  the 
partial  derivatives  of  Z  evaluated  at  the  steady-state 
operating    point. 
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B.   LINEARIZATION  OF  THE  MODEL  EQUATIONS 

There  are  basically  four  equations   to   linearize   for 
the  two-degree-of- freedom  model: 

ZDDOT  =  F  (Z,THEIA, MB) 

THETADDOT  =  G  (Z  ,  THET  A,  TH  ST  ADOT  ,  MB) 

M3D0T  =  H  (ZrMB) 

P3BAH  =  I(Z,MB) 


therefore  the  equatioas   for   the   differential   state 
variables  are: 


dZDDOT  =  (dF/dZ)  *dZ+  (d?/dTHETA)  *dTHETA+  (dF/dM3)  *dMB 

(III-1) 

iTKETADOT    =     (dG/dZ)  *dZ«-  (  iG/dTHETA)  *iTHET  A+  (dG/iM3>  *dMB 

(III-2) 

dM3D0T    =     (dH/dZ)  *dZ    +  (i:i/dTHETA)  *dTtfET  A+  (dH/dMB)  *dMB 

(III-3) 
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iP3BAR    =    (dl/dZ)  *dZ  +  (dl/dTHSIA)  *dIHETA  +  (dl/dMB)  *dMB 

(III-4) 

Note   that   the    left    hand    side    variables    are 
differentials,  not  time  derivatives. 

C.   LINEARIZED  EQOATIONS 


1.   Plenum  Air  Mass  FLdw  Derivatives 


Qin  =  N*(QiO-Ki*PBBA3) 


(11-12) 


Qout  =  Cn*Al* (2*PBBAR/RH3A)  **1/2 


(H-13) 


MBDOT    =    RHOA*(Qin-Qout> 

(11-14) 
as   a   result   the   differential   equation   for    3BD0r    becomes: 

dMBDOT    =    RH0A*(dQin-d2Ddt) 

(III-5) 

dQin  =  (dQin/iZ)  *dZ<- (dQin/dPBBAR)  *dPBBAR  +  ndQin/i^B)  *dMB 
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(III-6) 


dQin/dPBBAH  =  -N*Kq 

The  partial  derivatives  of  Qin  with  respect  to  all  other 
variables  are  zero. 


dQout  = 

(dQout/dZ)  *dZ  +  dQout/dPBBAR)  *dPBBAR+ (dQour/dMB)  *dMB 


dQout/dPBBAR  =  Cn*Al*  ( (3HOA/2*PBBAB)  **1  /2)  /RHO& 

(III-7) 

The  partial  derivatives  of  Qout  with  respect  to  all 
other  variables  are  zero.  Inserting  aquations  EII-6,and 
III-7  into  11-14  yields  the  equation  for  the 
differential   variable   dSBDOT: 

dMBDOI    =    (-RHOA*N*Kq-Ca*M*  (  (RHOA/2*  PBB  AR)  **1/2)  )  *dPBBAR 

(11-15) 

2  •       Plenum   Gau^e   Pressure   Dec  i/  at  i  ve.s^ 

Pb    =   Pa*(Hb/RHOA*Vb) **GAMHA 

(III-8) 
and 

PB6&R    =    Pb-Pa 
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Since  Pa  is  considered  constant  within  the  time  frame  of 
this  model, 

dPBBAR  =  dPb 

(III-9) 
is  a  result, 

dP3BAR  = 

GAdMA*Pa*  (MB (0) /RH0A*Vb (0)  )  **  (GAMMA-1)  *d (MB/RH0A*VB) 

where 

d  (MB/RHOa*VB)  =  (1/RH0A*V3  (0)  )  *dMB+ (MB/RHOA)  *d(1/VB) 

d(M3/RH0A*7b)  =  (1/RHO&*73  (0)  )  *dM3- (MB/RHOA*V3**2)  *dVB 

dV3  =  d(VN-AB*LD) 

dV3  =  -A3  (0)  *dLD/dZ-LD  (0)  *dA3/dZ 
remembering  that: 

LD  =  Z+ZS 
then 

dLD/dZ  =  1 

we  now  need  the  partial  of  effective   plenum   roof 
area  with  respect  to  z.   Defining 
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A3    =    ABtfL- (A3WL-ABR) *LD/3UBHGT 

(111-10) 

ABWL    =    HIDTH*LPML 

ABa    =    WIDTH*LPR 
farther   defining 

ABDIFF    =    ABWL-A3E 

which      is    a   constant,    t h s    effective    planum   roof    area    is: 


A3    =    ABWL-ABDIF?*LD/BQBH3r 

using    the    previous    definitions   the    partial   of      affective 
roof   area    becomes: 


dAB/dZ    =   -ABDIFF/BUBHGT 

(111-11) 
and   the   differential   voLume    is: 

d?b    =-  (AB  (O)-LD  (0)  *2*(A3DIFF/BUBKGT)  )  *dZ 

(111-12) 

with    the   result    that 

d  (.13/HHOA*Vb)     =     (1/RHOa*7b  (  0)  )  *dZ    + 
(Mb  (0)  /3HOA*Vb  <0)**2)*(AB(0)-2*LD(0)  *  (ABDIFF /BEJBHGT)  )  *dZ 


39 


finally , 

dPBBAR  =  GAMMA*Pb(0)  *((1/Hb(0)  )  *dHb  + 

(1/Vb(0)  *  (AB  (0)  -2*LD  (0)  *  (ABDIFF/3  0B3GT)  )  )  *dZ 


(111-13) 


3 •   Heave  Force  Derivatives 

The  total  heave  force  equation  is: 

H  (total)  =  HBF+HB&+HSF+HS&+HPLAN+HPBES 

The  total  partial  ierivatives  will  be  done  term  by 
term. 

For   the   forward   buoyant   force   the    following 
simplifications  are  made: 

HBF  =  -2*RH0*G*L1*LDBAR1* (LDB AR1 / (2* TAN (DR1)  )  +W510) 
let  K1  =  2*RH0*G*L1  then, 

HBF  =  -K1*LDBABT*(LDBAR1/(2*TAN  (DR1)  )  +WS10) 

where    LDBAR1    =    (LD-L5TAN  (TH2TA)  )  .    Combining   the 
previous  two  ejuations, 

HBF  =  -K1*  (LD-15*TAN  (THETA)  )  * 

(  (LD-L5*TAN  (THETA)/  (2*TAN(DR1))*tfS10) 
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assuming  small  pitch  angles,  TAN  (THETA)  can  be 
approximated  by  THETA  tfhich  resul-s  in  the  following 
simplification: 

HBF  = 

-K1*(LD-L5*THErA)  *  (  (LD-L5  *THETA)  /(2*TAN  (DR  1) )  +»5  10) 

Taking   the   partials   of   HBP    with,    respect   to      Z      and 
THETA    yields: 

d(HBF)/dZ    =    -K1*  ((LD  (0)  -L5*THETA(0)  /TA»  (DB1) +WS10 

(III- 14) 

and 

d  (HBF)  /dTHETA    = 

K1*L5*  ((LD{0)  -L 5* THETA  (0)  )  /  (2*TAN  (D/r1 )  )  ♦  WS10) 
+  K1*  (-L5/(2*TAN  (D31)  |  *(LD(0)  -L5*  THETA  (0)  ) 

Which    simplifies    to    : 

d  (HBF)  /dTHETA    = 

K1*L5*  (  (LD  (0)  -L5*THETA(0|  )  /(TAN  (DR1)  )  +H5  10) 

(111-15) 

Using   the    same   approach,    the    partials   for    buoyancy 
force    aft    are: 

3  (HBA)  /dZ    =    -K2*  ((LD  (0)  +L6*THETA  (0)  )  /(TAN  (D32)  )  «-WS20) 

(111-16) 
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d  (HBF) /dTHETA    = 

-L6*K2*  ((LD(O)  +L6*THET&  (0)  /(TAN  (DR2)  )  +WS10) 


(111-17) 


For  the  forward  seal  the  total  differentials   are: 

DHSF  =(d(HSF)  /dZ)  *dZ+(d  (HSF) /dTHETA)  *dTHET 
+  (d{HSF)  /dPBBAR)  *dPBBAR 

the  following  simplification  is  made; 

HSF  =  -P3BAR*WIDTH*  (LD-L3*TAN  (THSTA)  /SIN (31) 

(11-19) 
let  WIDTH/SIN  (31)  =  K3  and  TAN  (THETA)  =  THETA,  then, 

HSF  =  -K3*P3BAR*  (LD-L3*IHETA) 

with  the  resultant  partials  being: 


d(HSF)/dZ  =  -K3*PBBAR(0) 


(111-18) 


i  (HSF)  /dTHETA  =  K3*L3*PBB  AR  (0) 


(111-19) 


d  (HSF)  /dPBBAR  =  -K3*  (LD  (0)  -L3*THETA  ( 0)  ) 


(111-20) 


42 


The  partials  of  after  seal  force  are  founi  in   the 
same  manner  to  be: 

dHSA  =  (d  (HSA)/dZ)  *dZ+(d(HSA)  /dTHETA)  *dTH2TA 
+  {d(HSA/dPBBAR)  *dPB3AR 


HSA  =  -  (PBBAR+PDIFF)  *MIDTH*(LD+L4*TAN  (THSTA)/SIN  (32) 

(11-20) 
let  K4  =  WIDTH/SIN  (32)  and  assume  small  IHETA  so  that; 

HSA  =  -K4*  (PBBAR  +  PDIFF)  * (LD  +  L4*THETA) 
the  partials  then  becoae: 


d  (HSA)  /dZ  =  -KU*  (PBBAR(3) -t-PDIFF) 


(111-21) 


d  (HSA)  /dTHETA    =    -K4*  (PB&3  (0)  +PDIFF)  *L4- 


(111-22) 


a  (HSA)  /d?B3AR    =    -iU*  (LD  (0)  +L4  *THETA  (0)  ) 


(111-23) 


Next    the   planing   force    partials   are    derived: 


HPLAH    =    -2*PLCDEF*KHO*43EA*V**2*PI*THETA 


let    KPLAN    =     (2*?LC0EFF*RH0*AREA*PI)     then, 
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HPLAN  =  -KPLAN  *V**2*THETA 
and, 

d(HPLAN)  =  (d  (HPLAN)  /dlHET  A)  *dTHETA 

for  this  constant  speed  model,  therefore, 

d  (HPLAN)  /dTHETA  =  -KPLAN  *V  (0)  **2 

(111-24) 
Plenum  Gauge  Pressure  lift  was  defined  to  be: 

HPRS5  =  -AB*PBBAR 

(II-11c) 
as  a  result: 

a  (HPRES) /dPBBAR  =  -A3(0) 

(111-25) 

4 •   Pitch  Moment  Partials 

The  total  moments  acting  on  the  craft  are: 

PTOTAL  =  PBF  +  P3A  +  PSF+P5A*PPLAN  +  PPR2S  +  ?DA!1P 

(11-38) 

The    partials    of    11-38    will    be    found    term    oy      term, 
Deginning    with    the    forward   buoyant    moment: 
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P3F  =  -HBF*L5 

(11-23) 

therefore: 


d  (PBF)  =  -L5*d(HBF) 
so  that: 

d(PBF)  =  -L5*  (d  (HBF)  /dZ)  *dZ-L5*(  (d  (HBF)  /dTHETA)  *dIHETA 

(111-26) 

and  the  partial  of  the  after  buoyant  moment  is: 


PBA  =  L6*HBA 


(11-26) 


dP3A   =    L6*dHBA 
therefore , 

d(?BA)  =  L6*  (d  (HBA)  /dZ)  *dZ+LS*  (d  (HBA)  /dTHETA)  *dTHETA 

(111-27) 

Seal  forces  partials  must  include  the  dynamic 
variation  in  lever  arm  length,  which  for  tha  forward 
seal  is: 

PSF  =  -HSF*PL5F 
or, 
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PSF    =     (PBBAR    *ASEAL1) *PLSF 
therefore , 

d  (PSF)  /dPBBAR    =    ASSAL1  (0)  *PLSF  (0) 

and 

d  (PSF)  dZ  =  PBBAR  (0)  *d  (PLSF*ASEAL  1)  /dZ 
where 

d(AS2AL1*PLSF)  /dZ  =  PLSF  (0) *d  (ASEAL1 ) /dZ  + 
ASBAL1  (0)  *d(PLSF)/iZ 

ASEAL1  =  WIDTH*  (LD-L3*T SET  A) /SIN  (31) 
and 

PLSF  =  L3-(LD-L3*IHETA)  /2*SIN(31) 
therefore , 

i(?LSF)/dZ  =  -1/(2*SIN  (31) 

and 

i(ASSAL1)/dZ    =    KIDTH/SIN  (3  1 ) 
as    a   result, 

d  (PSF)  /dZ    =    PBBAR  (0)  *(  (PLSF(O)  *WIDTH) /SIN  (31; 
ASEAL1  (0)  /{2*SIN  (31)  ) 


(111-23) 
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(111-29) 


a  (psf)  /dTHETA    =    PBAR  (0|  *a  (PLSF*ASEAL1)  /dTHETA 


d  (PLSF*ASEAL1) /dTHETA  = 

PLSF  (0)  *d  (ASEAL1)/drHETA+ASSAL1  (0)  *d  (PLSF)  /dTHETA 


d  (ASEAL1)  /dTHEIA  =  -L3*»IDTH/SIN  (31 ) 

d  (PLSF)  /dTHETA  =  L3/ (2*SIN  (3  1 ) 
as  a  result, 

d  (PSF)  /dTHETA  =  -PB3&R  (0)  *  (  (PLSF  (0)  *L3*  WIDTH) /5  IN  (3  1  )- 

ASEAL1  (0)  *L3/(2*SIN  (31)  ) 

(111-30) 

For    the   after    seal    the   derivatives    proceel    in      the 
manner: 

PSA    =    -{PB3AR    +    PDIFF)  *ASEAL2*PLSA 

d  (PSA) /dPBBAR    =    -ASEAL2  (0)  *PLSA(0) 

(111-31) 

In  the  same  manner  as  shown  in  the  previous 
paragraph,  the  partials  of  moment  iua  to  seals  aft  with 
respect   to   the    other    variables   are: 
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d  (psa)  /dZ  =  -PBBAR  (0)  *  (  (  PLSA  (  0)  *HIDTH)  /SIN(32)  + 
ASEAL2  (0)  /(2*SIN  (32)  ) 


(111-32) 


d  (PSA)  /dTHETA    =    -PBAR  (0)  *  (  (PLS  A  (0)  *L4*WIDTH)  /SIS)  32)  + 
ASEAL2(0)  *L4/(2*SIN  (32)  ) 

(111-33) 

The  partial  of  planing  moment  is: 


PPLAN  =  L7*HPLAN 


(11-30) 


d  (PPLAN)  =  L7*d(HPLAN) 


d  (PPLANJ/dTHETA    =    -L7*KPL AN*V (0)  **2 

(111-34) 

For    this   constant   speed   model   all    other    partials  are 

identically      zero.         Th=      Plenum   Gauge   Pressure    also  has 

dynamic   variation      in      the      moment      lever      arm      and  the 
derivatives   are: 


PPRE3    =    -XCP*HPRES 

(11-34) 

therefore , 

dPPRES    =    -HPRES  (0) d  (XCP)  -XCP  (0) *d  (HPRES) 
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where 

XCP    =    XCPO+  (LDIFF*LD)/(2*BUBHGT) 
therefore, 

d(XCP)/dZ    =    LDIFF/(2*BUBHGT) 
and 

d(PPRES)     =    -  ({HPRES  (0)  *LDIFF)  /  (2*BU3HGT)  )  *dZ 
-XCP(O) *d(HPRES) 

substituting   equation    (111-25)    yields: 

d(PPRES)     =    -{  (HPRES(O)  *LDIFF)  /(2*BOBHGI)  )  *dZ 
+XCP{0) *AB (3) *dPBBAR 

Finally,    the   partial   of    pitch   damping   moment   is: 

PDAHP    =    -2*PI*XLSS**2*YLSS**2*RH3*u*THETADOT/8 

(H-37) 

let 

2*PI*XLSS**2*YSS**2*RHQ/8    =    KDAMP 

then 

PDAtfP    =    -KDAMP*V*THETADDT 

as    a   result,    since    this    is    a    constant    speed   modal, 
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a  (PDAKP)/dTHETDT  =  -KDA!1P*V(0) 

(111-35) 

where   THETDT   is   the   first   derivative   of  theta  with 
respect  to  time. 


5.   Total  Dif f ergntials 


Since  F  =  M  *  A  and  the  interest  is  in 
acceleration  terms,  A  =  F/.1  and  dA  =  d(F)/M.  the  result 
being  that  all  differentials  mast  oa  divided  by  craft 
mass  to  achieve  the  lesired  linearized  heave  results. 
The  same  general  coaanent  holds  true  for  moment 
differentials,  that  is;  they  all  must  be  divilsd  by  the 
craft  moment  of  inertia  about  the  Y-axis. 

All  dif f arentials  must  be  sorted  by  the  variable 
of  differentiation  so  that  the  linearized  differential 
equations  of  motion  may  be  written.   As  an  example; 

dZDDOT  =  ((d(HBF)/dZ)  *■  (  d  (HBA)  /dZ)  )  /MASS)  *dZ  + 

The  State  Variable  representation  of   the   system 

will   be   presented   in  a  later  section  of  this  chapter, 

however  it  is  convenient  to  define  the  State   Variables 
at  this  time. 

t 

i  =  (dzrazDoi,aiHErA,dra£TDT,dHB) 

(rii-36) 

and 
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Y  =  dPBBAR 

(111-37) 

since   PBBAR   is  not  a  state  variable  the  definition  for 
dPBBAH  must  be  substituted  for  dP33AR  in  all   equations. 

6 •   Total  Heave  Force  D iff erentials  in  Z 

Here  the  total  differentials  will  be  found.  No 
explanatory  text  will  be  inclided  unless  special 
suostitut ions  requiring  an  aside  are  necessary. 

The  buoyancy  terms  are: 

DH3Z  =  d(HBF)/dZ  +  d(HB.\)/dZ 

(111-38) 

A  term  must  be  included  which  was  not  previously 
derived  to  account  for  the  sensitivity  of  seal  forces 
due  to  PBBAR  since  it  is  not  a  State  Variable. 

DHSPZ  =  (d  (HSF)  /dPBBAH)  *  (d  (PBBAR)  /dZ) 

as  a  result, 

d  (HSF)  /dPBBAfi  =  -K3*  (LD  (0)  -L3*THETA  (0)  ) 


(111-20) 


and 


DHSFPZ  =  -K3*  (LD  (0)  -L3*THETA(0)  )  *d  (PBBAR) /dZ 
similarly  for  the  aft-ec  seal  moment: 
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d  (HSA)  /dPBBAR  =  -K4*(LD(D)  +L4*THETA  (0)  ) 

(111-23) 
aad 

DH3APZ  =  -K4*(LD  (0)  +L4*rHErA(0)  )  *d(PBBAH) /dZ 

The  total  seal  force  differential  in  Z  becomes: 

DHSZ  =  d(HSF)/dZ  ♦  d(HSA)/dZ  +DHSFPZ  +  DKSAPZ 

(111-39) 


For  plenum  gauge  pressure   lift   differential   the 
same  substitution  must  be  made. 


dHPRES  =  -AB(0)  *dPBBAR 
therefore, 

DHPZ  =  -AB  (0)  *d  (PBBAR)  /dZ 

(111-40) 

Define  the   sensitivity  coefficient  of  ZDDOT  to  Z 
to  be: 

DZZ  =  (DHBZ  +  03SZ  +■  DHPZ) /MASS 

(111-41) 

7.   Total  Differentials  in  THETA 
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The    following    partials    are   now    defined: 


DH3TH    =    d  (HBF)  /iTHETA    +    d (HBA) /dTHETA 


(111-42) 


DHSTH  =  d  (HSF)  /dTHETA  +  i  (H3A)  /dTHETA 


(111-43) 


DH?TH  =  d  (HPLAN) /dTHETA 


(111-44) 


Define   the    sensitivity    coefficient    of    ZDDDT    due   to 
IHETA    to    be; 


DZTH    =     (DHBTH    +■    DHSTH    +    DHPTH)  /MASS 


(111-45) 


3-   Total  Differentials  due  to  M3 

Here  since  dPB3A3  is  not  a  state  variable  the 
partial  of  PBBAR  with  respect  to  MB  must  be  suostituted 
in  the  plenum  lift  and  seals  lift  terms. 

3  (HSF)  /dPBBAR  =  -rK3*  (LD  (0)  -L3*THETA  ( 0)  ) 

(111-20) 
by  substitutiug  d(PB3AR|/dMB  for  dPBBAR  the  result  is: 
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d(HSFP)/dMB  =  -K3*(LD(3|  -L3*rHETA(0)  )  *d  (PBBAR) /dM3 

(111-46) 
similarly, 

d(HSAP)/dMB   =    -K4*  (LD  (0)  +L4*THETA  (0|  )  *d  (PB3AR)  /dMB 

(111-47) 

d(HPSES)     =    -AB(O)  *dPBBAR 

using  the  previous  substitution  yields: 

d(HPBES)/dMB  =  -AB (0)*i (PBBAR) /dMB 

(111-48) 

The  sensitivity  coefficient  of  h^ave  due  to  MB  can 
now  be  defined  to  be: 

DZMB  =  (d  (HSFP)  /dMB  +-d  (35  AP)  /dMB  +  1  (HPBES)  /dMB)  /MASS 

(111-49) 

9 .   Total  Pitch  Mo  sent  Differentials  in  Z 

DP3Z  =  d(PBF)/dZ  +  d(PB4)/dZ 

(III -50) 

The   dPBBAR   substitution   yields   the   partial  of 
seals  with  respect  to  pressure  with  respect  to  Z. 
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DP5FPZ  =  ASEAL1  (0)  *PLSF  (0)  *d  (PBBAS)  /dZ 


(111-51) 


DPSAPZ  =  -ASEAL2(0)  *PLS&(0)*d  (PBBAR)  /dZ 

(111-52) 

therefore  the  total  seal  differential  in  Z  becoaes: 


DPSZ  =  DPSFPZ  +  DPSAPZ  +  d(PSF)/dZ  +  d(?SA)/dZ 


(111-53) 


Next    the    Plenum   3a.ige   pressure    term    is   isrived: 
d(PPRES)/dZ    =    -HPRES  (0)  *LDIFF/(2*B0BHGr) 

d  (?PRES)/dPBBAR    =    AB  (0)  *XC?  (0) 
therefore, 

DPPRESZ   =    AB  (0)  *XCP(0)  *i  (PBBAR)  /dZ 
and 

DPRESZ    =    d(PPP.ES)/dZ    +  d(  PPRES  Z)  /dZ 

(111-54) 

The      sensitivity    coefficient    of   IHETA   with    respect 
to    Z    becomes: 

DTHZ    =     (DPSZ    ♦    DPSZ     *    DPPRES)/LYY 
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(111-55) 

10.   Total  Moment  Differentials  in  THETA 

DPBTH  =  (d  (PBF)  /dTHETA  +  d  (PBA  )  /iTH3  TA) 

(111-56) 

DPSTH    =    (d  (PSF)  /dTHETA    <•    d  (PS A)  /dTHETA 

(111-57) 

DPPLTH  =  d(PPLAN) /dTHETA 

(111-58) 

define  the  sensitivity  coefficient  for  THETA  to  be: 

DTHTH  =  (DPBTH  +  DPSTH  *  DPPLTH)  /HI 

(111-59) 

Also,  define  the  sensitivity  coefficient  in  THETDT 
to  be : 

DTHDTH  =  (d  (PDAMP)  /dTHDT)  /IYY 

(III -60) 

71.       Total    Moment    Differentials    in    MB 
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dPSAMB    =     (d  (PSA)  /dPBBAR)  *  (d  (PBBAH)  /dMB) 

(111-61) 

dPSFMB    =(d  (PSF)  /dPBBAR)  *  (d  (PBBAR)  /dMB) 

(111-62) 

DPSMB  =  DPSAMB  +  DPSFMB 

(111-63) 

Inserting  the  dPBBAH  substitution  yields  the  partial   of 
plenum  gauge  pressure  moment  with  respect  to  M3: 

d  (PPSES)/dPBBAR  =  -XCP  (0)  *DHPR3S 

d(HP2ES)/dMB  =  -  AB  (0)  *i  (  PB3AR)  /dMB 

DPPMB  =  XCP(O)  *AB(0)  *d(?33AR) /dMB 

(111-64) 

Define  the  sensitivity  coefficient  of  THETA.  due  to 
M3  to  be: 

DTHMB  =  (DPSMB  +  DPPMB)  ZIYI 

(111-65) 

12.       Total   Differentials    of    Mr    Mass 
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d(tfB)/dM3  =  (d  (MB) /dPBBAR)  *  (d  (PBBAR)  /dMB) 
define  the  sensitivity  coefficient  to  be: 

DMBMB  =  d(MB)  /dKB 

(111-66) 
define  the  sensitivity  coefficient  in  Z  to  be: 

DM3Z  =  d(MB)  /dZ 

(111-67) 


D.   TOTAL  MODEL  SENSITIVITY  COEFFICIENTS 


1  •   Heave  Sensitivity  Coefficients 

The  sensitivity  coefficients   for   the   linsarized 
differential  equation  in  heave  acceleration  are: 

DZZ,  DZTH,  and  DZHB. 

2 •       Pitch   Sensitivity.  Coefficients 

The    sensitivity   coefficients      for      the      Linearized 
differential   equation   in    pitch   acceleration   ara: 

DTHZ,    DTHTS,    DTHDTH,     and    DTHHB. 
3.       Plenum    ^ir    Mass   Sensitivity    x2§££ isi^Si§ 
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The  sensitivity  coefficients  for  plenum   air   mass 


are: 


DMBZ,  and  DMBMB. 

4 •   Plenum  Gauge  Pressure  Sensitivity  Coefficients 

.  The   sensitivity   roefficients   for  PL  a  nun   Gauge 
Pressure  are: 

DPBZ,  and  DPBMB. 

E.   STATS  SPACE  REPRESENTATION 

The  state  matrix  form  of  a  dynamic  system  is: 

XDOT   =  AX  +  BU 

where  X  =  a  vector  of  variables  representing  the  states 
of  a  set  of  first  order  linear  differential  equations, 
and  0  is  fc^e  vector  of  forcing  functions,  and  A  and  B 
are  the  matrices  of  the  respective  gains.  The  form  of 
the   equations   is: 

X1D0T  =  X2 
X2D0T  =  X3 
X3DOT    =    K1*Xt    +    K2*X2    *    K3*X3 
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for      the      linearized      model      of      the      CAB      SSS    the   state 
variables    were   defined   to    be    : 


X    =    (dZ,    dZDOT,    dTHETA,    dTHETDT,    daB) 


As      a      result      the      state   matrix    representation    of    the 
linearized    unforced   Heave-Pitch   model   is: 


X1DOT 


X2DOT 


X3D0T 


X4UOT 


X5D0T 


DZZ 


DTHZ 


DMBZ 


DZTH 


DTHTH         DTHDTH    DTHM3 


0  "" 

Xl" 

DZHB 

X2 

0 

X3 

DTHMB 

xu 

DMBMB 

X5 

and 


dPBBAR  =  DPBZ*X1  +  DPBMB*X5 
111-69 


(111-68) 
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IV.   JIME  fiOMAIN  VALIDATION 


The  linearized  model  developed  in  chapter  III  was 
validated  against  the  the  Six-degree-of-freedom 
nonlinear  model.  There  were  several  changes  made  in  the 
simulation  techniques  used,  and-  at  the  same  time  the 
changes  in  center  of  pressure  dynamic  variation  and 
planing  forces  were  entered  into  the  6  DDF  model. 
Establishment  of  eguibibrium  for  the  linear  model  was 
accomplished  off-line  using  a  TI-59  programmable 
hand-held  calculator  with  printer. 

The  first  change  in  simulation  technique  is  the 
operating  point  at  which  the  linear  sensitivity 
coefficients  were  evaluated.  In  Mclntyre's  model  these 
coefficients  were  evaluated  at  the  initial  operating 
point.  In  this  work,  the  coefficients  are  evaluated  at 
the  final  steady-state  equilibrium  point  since  this  is 
customary  in  linearization  work. 

The  second  major  change  in  simulation  was  the  method 
usad  to  determine  the  planing  force  loss  coefficient 
(PLCOEFF)  and  the  establishment  of  the  arbitrary 
location  of  the  initial  plenum  center  of  pressure.  In 
this  model  all  heave  forces  are  well  defined  (as  nodeled 
here)  with  the  exception  of  planing  force  lift.  As  a 
result,  craft  weight  and  all  heave  forces  axcapt  that 
due  to  planing  lift  are  summed  and  the  residual  assumed 
to  oe  due  to  the  planing  phemonema.  This  residual  is 
then  used  to  determine  (off-line)  the  value  of  PLCOEFF. 
Note  that  this  method  also  accounts  for  any  lift  due  to 
aerodynamic   forces   on   the   bow  that  were  disregarded. 
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Next  all  moments  are  sunaed  about  the  C.G.  except  the 
moment  due  to  plenum  pressure  lift.  The  residual  moment 
is  assumed  to  be  due  to  plenum  pressure.  This  residual 
is  used  (again  off-line!  to  establish  the  initial  center 
of  pressure.  This  accounts  for  any  moments  not 
previously  accounted  for  (such  as  aerodynamic  forces  on 
the  bow.)  Mclntyre  used  the  planing  lift  coefficient  to 
force  the  model  into  agreement  with  the  6  DOF  uodel  for 
natural  pitch  frequency  and  used  a  large,  arbitrary,  and 
static  moment  to  force  the  models  into  steady  state 
pitch  angle  agreement.  It  was  felt  that  using  the 
present  method  of  determination  would  laave  the 
equations  in  a  form  that  would  be  more  amenable  to 
extraction  of  the  factors  most  directly  affecting  craft 
characteristics  as  modeled  herein. 


LINEAR  PROGRAM  DESCRIPTION 


User's  instructions  have  been  included  as  appendix  A. 
The  program  flow  can  be  isolated  into  four  main 
sections. 

1.  First  the  values  for  craft  weight,  draft,  pitch 
angle,  plenum  pressure  and  speed  are  entered.  To 
initialize  plenum  air  mass,  draft  is  used  to  solve  for 
initial  plenum  air  volume,  which  is  then  used  with 
initial  plenum  gauge  pressure  to  solve  equation  IT-15 
for  plenum  air  mass.  Sguation  11-12  is  then  solved  for 
Sin.  Qout  is  set  equal  to  Qin  and  equation  11-13  is 
solved  for  seal  leakage  area.  Ail  other  heave  and  pitch 
nonlinear  equations  are  solved  and  their  values  are 
listed  along  with  residuals  of  heave  aid  pitch 
summations.  These  resiiuals  are  added  to  the  initial 
values  of  C.G.   and  angalar  acceleration  and  tne  initial 
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values  are  saved  for  the  time  domain  simulation  section. 

2.  The  same  operations  are  carried  out  for  tie  steady 
state  values  of  weight,  draft,  pitch  angle,  and  plenum 
pressure.  The  residuals,  forces  and  moments  are  again 
listed. 

3.  These  steady  state  values  are  used  to  evaluate 
the  sensitivity  coefficients  developed  in  chapter  III 
which  are  also  listed. 

4.  Finally  the  stats  variable  solution  is  simulated 
using  a  four-point  fixed  step  Runge-Kutta  algorithm  and 
listed  for  the  print  interval  desired.  The  option  is 
included  to  output  either  printer  plot  grapnicaL  display 
using  the  NPS  printer  plot  suoroutine  PLOTP  or  the  X-Y 
Versatec  plot  subroutine  OHAWP.  The  latter  should  be 
used  for  final  runs  only,  with  PLOT?  usad  on 
intermediate  runs.  Provisions  are  made  for  up  to  300 
discrete  output  points  for  the  time  domain  tabular 
listing . 


3.   SIX-DEGREE-OF-FREEDOM  MODEL  SOLUTION 


Two  changes  were  made  in  the  6  DOF  model  to  account  in 
a  more  analytic  manner  for  the  observed  tendancy  of  the 
XR-3  to  pitch  down  as  speed  increased. 

The  first  change  was  to  remove  an  empirical  eguation 
used  to  correct  the  center  of  pressure  L5]  and  replace 
it  with  the  dynamic  center  of  pressure  variation 
developed  in  chapter  II. 

The  second  change  rfas  to  include  in  the  sidewall 
subroutine  the  nonlinear  equation  for  planing  forces, 
with   the   loss   coefficient   adjusted   to   make   up  (as 
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closely  as  possible)  foe  lift  loss  due   to   the   revised 
modeling  of  the  plenum  chamber. 

Finally,  the  6  DOF  modal  was  ran  at  both  6722  lbs  and 
6050  lbs  weight  to  establish  the  initial  and  final 
operating  points,  followed  by  the  weight  removal 
transient  run  for  time  domain  validation  of  tha  linear 
heave-pitch  model. 

It  should  be  noted  that  furtnar  study  of  the  6  DOF 
model  is  needed  to  evaluate  tha  results  of  the  two 
changes  and  to  more  closely  establish  the  arbitrary 
initial  location  of  the  center  of  pressure.  The  large 
amount  of  CPU  and  turn-around  program  times  raguired  to 
run  the  nonlinear  model  precluded  including  that  study 
in  this  work.  Also,  for  tiie  nonlinear  model, 
X-direction  thrust  was  held  constant  during  the  run  with 
the  result  that  speed  was  30.33  knots  and  still 
increasing  at  the  end  of  the  five  second  rua.  This 
effect  will  be  further  liscussed  in  tha  next  section. 


C.   COMPARISON  OF  RESULTS 

The  graphical  output  for  the  five  second  runs  are 
shown  in  figures  5  through  3  tor  the  linear  aodel  and 
figures  9  through  12  foe  the  nonlinear  model.  The  fixed 
integration  step  size  was  0.005  se:  for  the  linaar  model 
and  0.001  sec  for  the  nonlinear  modal. 

1 .   Planum  Gaugja  Pressure  Transient 

(see  fig.  5  and  9| 
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In  general  both  models  exhibited  a  rapid  initial 
reduction  in  pressure  fallowed  by  a  smooth  transient 
back  towards  the  final  operating  pressure.  For  the 
linear  model  the  minimum  occur  ad  at  100  milli-seconds, 
and  110  milli-seconds  for  tha  5  D3F  model,  with  the 
respective  pressures  being  21.38  and  21.30.  Tie  final 
values  were  24.31  and  24.08  for  the  linear  and  nonlinear 
model   respectively. 

As  can  be  seen  from  the  graphical  data  excellent 
agreement    was    acheived   between   the    two    models. 

2 .  Draft   Transient 

{see  fig.  6  and  13) 

Both  models  exhibit  a  smooth  exponential  type 
transient  towards  the  final,  lower  steady-stata  draft. 
Tha  linear  model  transient  settles  towards  a  lower  value 
(6.04  inches  as  compared  to  6.44)  which  is  in  part 
attributed  to  a  difference  in  modeling  of  the  seal 
forces.  The  nonlinear  aodel  uses  a  complex  combination 
of  hydrodynamic  and  hydrostatic  lift  where  the  linear 
model  uses  a  simple  pressure  differantial  effact.  Mso, 
since  the  nonlinear  modal  is  settling  to  a  significantly 
lower  pitch  angle  (see  para.  4  below) ,  the  planing  force 
effect  is  reduced  in  that  model. 

3 .  C .G  Acceleration 

(see  fig.  7  and  11) 

Both  modals  show  a  rapid  positive  transient  in 
acceleration  that  closely  follows  the  negative  transient 
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in  plenum  pressure.  Both  models  reached  the  initial 
peak  positive  C.G.  acceleration  at  approximately  20 
milli-seconds  after  the  plenum  pressure  completed  the 
negative  transient.  The  linear  model  appears  to  have 
greater  coupling  between  pitch  and  heave  looking  at  the 
C.G.  acceleration,  but  this  is  not  actually  supportable 
since  the  linear  model  has  greater  pitch  amplitude 
excursions  about  the  exponential  pitch  transient  to  the 
lower  steady  state  value  (see  fig.  8  and  12.)  In 
general  the  two  models  show  excellent  agreement  in 
transient  and  damped  C.  3.  acceleration  behavior. 

4  •   Pitch  Tr a^nsien t 

(see  fig.  8  and  12) 

It  is  in  this  area  that  the  two  models  are  in 
greatest  disagreement.  Both  exhibit  a  general  damped 
sinusoidal  behavior  about  a  center  exponential  transient 
to  a  lower  steady  state  value.  In  both  :isas  the 
initial  excursion  is  downward.  The  natural  pitch 
fequency  of  linear  model  is  approximately  4.3  radians 
per  second,  whereas  the  nonlinear  model  exhibits  a 
natural  pitch  frequency  of  6.98  radians  per  second. 
There  is  felt  to  be  one  major  factor  contributing  to 
this  large  discrepancy,  which  is  the  inclusion  of  the 
planing  force  in  the  nonlinear  model.  Data  presented  by 
Thaler  and  GerbaO'G  showed  a  natural  pitch  frequency  of 
5.3  radians  per  second  and  much  higher  oscillations 
about  the  exponential  transient.  It  is  felt  that  the 
simplified  planing  forces  of  the  linear  nodel  are 
included  in  the  much  more  complex  hydrodynamic  equations 
of  the  nonlinear  model  and  are  therefore  redundant  in 
that  model.  Due  to  time  constraints,  that  avenue  of 
investigation   is   left   for  future  studies.   The  linear 
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model  is  tending  towards  a  steady  state  value  on  the 
order  of  0.45  degrees  whereas  the  nonlinear  model  is 
tending  towards  approximately  0.30  degrees.  This  large 
discrepancy  is  partiaLLy  due  to  the  redundant  planing 
force  in  the  nonlinear  model  and  the  fact  that  speed  is 
also  increasing  in  the  nonlinear  model  which  tends  to 
accentuate  the  planing  force  and  force  the  craft  to 
pitch  down. 

In  the  next  chapter  a  signal  flow  graph  will  be 
developed  and  Mason's  3ain  Rule  applied  to  explain  the 
linear  model  characteristic  response  in  a  more 
analytical  fashion. 


67 


SIGNAL  FLOW  GRAPH  AND  MASON'S  GAIN  RULE 


A.   INTRODUCTION 


The  State  Variables  defined  in  Chapt  II  are  used  to 
develop  a  signal  flow  graph  of  the  craft  h=ive-pitch 
dynamics.  Mason's  Gain  rule  is  applied  to  the  signal 
flow  graph  to  derive  analytical  S-domain  polynomials 
expressing  the  transfer  functions  for  pitch  moment  into 
pitch  angle,  heave  acceleration  into  draft  and  the  two 
cross  gains.  Using  typical  gains  calculated  in  the  time 
domain  program  developed  in  Chapt  III  the  roots  of  these 
polynomials  are  located  and  it  is  shown  that  pitch 
dynamics  may  be  very  closely  approximatsd  by  a 
second-order  system.  Heave  dynamics  are  also 
simplified. 

Using  a  simplified  approach  to  the  frequency  response 
of  the  craft  to  non  calm-water  conditions,  the  basics 
are  developed  for  later  sea-state  model  work.  A  dynamic 
system  modeling  program  available  at  tae  Naval 
Postgraduate  School  is  used  to  generate  Bode  plots  for 
the  craft's  frequency  response  in  Heave  and  Pitch  under 
the  simplifying  assumptions. 


SIGNAL  FLOW  GRAPH  ANALYSIS 


Fig     13    was    generated    using   the   sensitivity 


68 


coefficients  defined  in  Chapter  III.  From  that  graph  the 
following  feed-forward  and  feedback  loops  and  paths  were 
identified: 

1.    Feed   Forward,    U(THETA)     into    THETA. 

P1    =     (1/IYY)*  (VS**2) 

2.    Feed   Forward,    U(Z)     into   Z. 

P2    =     (1/MASS)  *  (1/S**2) 


3.     Feed    Forward,    0(T3ETA)     into    Z. 


P3    a     (1/IYY)  *DZIH*  (1/S**4) 


4.  Feed  Forward,  U(Z)  into  THETA. 


P4  =  (1/:iASS)  *DTHZ*(1/S**4) 


5.Feedbac)c    Loops. 


L1    =    DTHTH/S**2 


L2    =    DTHDTH/S 


L3    =     (DZIH)  *(DTHZ)  /S**% 


L4    =     (DZMB)  *(DMBZ)  /S  **3 
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L5    =    DZZ/S**2 


L6    =    DHBMB/S 


L7    =    (DZTH)*(DSBZ)  * (DT98B) /S**5 


C.       MASON'S    GAIN    ROLE    APPLIED 


Mason's    gain    rale    for    a    closed    multi-loop  system    is: 


T    =    SUM  (P)c*DELIAk) /DELIA 

whara , 

T   =   transmission    gain    from      an      input      to      in      output 
node. 

DELTA   =    the   graph   datarminant 

DELTAJc    =   cofactor   of   the   !cth   path 

Pic   =   gain   of   the   diract    path   from   input   to   output 

Using      Mason's    gain   rula   the    graph    daterminant     (DELTA) 
is   found    to    be    of   the  form: 


D(s)     =     (35*S**5    +    B4*s**4    +    33*s**3    +    B2*s**2    +    B1*s    + 
B0)/S**5 

The  forward  direct  path  0  (THETA)  to  THETA  is: 


N1(s)  =  (1/IYY)  *  (1/S**5)  *  (A3*s**3  +  A2*s**2  +  A1*S  +  AO) 


The  forward  path  0  (Z)  co  Z  was  found  to  be: 


70 


N2(S)  =  (1/MASS)  *  (1/S**5)  *(A3*S**3  *    A2*s**2  *  M*S  + 
AO) 

The  forward  direct  path  U(THETA)  to  Z  is: 

N3(S)  =  (1/IYY)  *DZIH*(S  +  A0)/S**4 

The  forward  path  for  draft  coupling  into   pitch   (U(Z) 
to  TH2TA  is: 

N4  =  (1/MA5S)  *DTHZ*  (S  +  A0)/S**4 


D.       SOLUTION    OF    THE    INDIVIDUAL    POLISOMIALS 


The    gains   used    below    were  taken    from      a      typical      time 
domain    run    as   developed    in   chapter   IV. 

DZZ    =   -2178.0 
DZTH   =    -121.3 
DZtfB    =    -5499. 0 
DTHZ    =    28.78 
DTHTH    =    -22.  3  1 
DTHM3    =    72.1 
DMBZ    =    -23.79 
DMBHB    =    -73. 1 1 

1 •       Solution    of    the   3ca£h    Determinant 

Inserting      the      nuaerical      values      for   the    typical 
sensitivity   coefficients    listed   above   and    using    only   the 
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dominant      terms      where      possible,    the    roots  of    the    graph 
determinant   were    found    to    be: 


DELTA    =     (S**2    +    71.8*S    f    2066)*(S    +    0. 5  086)  *  (S*  *  2    + 
0.72*S    +22.74} /S**5 

These  factors  represent  a  pair  of  compLex  poles 
located  at  -35.71+  J27.84  ,  a  real  pole  located  at 
-0.5085  ,  a  pair  of  -Duplex  conjugate  roots  located  at 
-0.360  +  j  4.755  ,  and  a  fifth  order  zero  located  at  the 
origin . 

2-      Solution   of   H-LIL 


N1  (S)     =    (1/IIY)*{S    +    0.427)*(S**2    *■    72.68*S    +    2146)/S**5 

This  polynomial  factors  into  a  real  zero  located 
at  -0.427  ,  a  pair  of  complex  conjugate  zeros  located  at 
-36.34  +  j  28.74  ,  and  a  fifth  order  pole  locate!  at  the 
origin. 

3-   Solution  of  M2.iSi_ 


N2(S)  =  (1/MASS)  *{5**2  f  0.7721*S  +  22.35)*(S  + 
72.38)  /S**5 

Which  factors  into  a  pair  of  complex  conjugate 
zeros  located  at  -0.3856  +  j  4.712  ,  a  real  zero  at 
-72.38  ,  and  a  fifth  order  pole  located  at  the  origin. 
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**■   Solution  of  N3iSL 

N3  (S)  =  (1/IYY)  *DZTH*(S  +  73.11)/S**4 

Which  represents  a  real   zero   at   -73.11,   and  a 
fourth  order  pole  at  the  origin. 

5-   Solution  of  N4JS^ 

N4 (S)  =  (1/MASS) *DIHZ*(S  +  73.11)/S**4 

Which   has   the  same   pole   and  zero  locations  as 
N3  (S)  . 

E.   SIMPLIFICATION  OF  TRANSFER  FUNCTIONS 


1 •  Sii il 

T1  (S)  =  THETA(S)  /0  (THETArS^) 

Therefore   from   Sason1 s   Gain   Rula: 

T1  (S)     =    »1  (S) /DELTA 


11  (SJ     =     (1/IYY)  *  (S    *    35.34    +    428.7)*(S    +    0.427)/ 

(S    +    35.1    +    j    27.84)*(S    +    0.503)*(S    *    0.36    +    j    4.755) 
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To  simplify  let  tie  complex  conjugate  zaros  at 
-36.34  +  J28.7  cancel  the  poles  at  -35.9  +  J27.84  ,  and 
the  real  zero  at  -0.427  cancell  the  pole  at  -0.508. 
The   result    is    the    simplified    second   order   system: 

T1  (S)     =    (1/IYY)  /  (S**2    ♦    0.72*s    +    22. 7U) 

Inspection  of  the  time  domain  pitch  response  with 
respect  to  natural  frequeny  and  damping  confirms  this 
basic    underdamped   second    order   response. 

2.       T2  (S) 


I2(S)     =    Z  (S)/0  (Z,5) 


I2(S)     =    (1/MASS)*(S    +    0.3856    +    j4.712)*(S    +    72.38)/ 

(S**2    +35.1    +    j    27.34)*(s    +    0.5035)*((S    +0.35    + 
J4.755) 

Canceling  the  pole  at  -0.36  +  j  4.755  with  the 
zero  at  -0.3856  +  J4.712  results  in  the  simplified 
response: 

P2(S)     =    (1/HASS)*(S    +    72.33)/(S    +    0.503)*(S**2    f    71.8*S 
+    2066) 

Further  simplification  can  be  made  if  (S  + 
72.33) /(S  +  0.508)  is  approximated  by  72.33/(3  +  0.508) 
(which  is  valid  for  longer  time  period  considerations.) 
As  a  result  T2(S)  becomes: 
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T2  (S)     =     (72.33/MASS)/(5*0.508)  *(S**2    *    71.8*S    +    2066) 

which  has  poles  located  at  -0.508,  and  -35.9+  J27.83 
This  is  in  quite  close  agreement  with  the  results 
obtained  by  Thaler  and  Serba  14  for  the  same  airflow 
condition  (Kq  =  1.0)  for  their  heave-only  linearized 
model.  Their  results  ware  a  real  pole  at  -0.33  and  a 
complex   pair    at    -35.1*    J34.4    .    This   is    of   the    form: 

T(S)     =    K/(S    +    0.38)*{S**2    +    70.2*S    +    2415) 


T3  (S)     =    Z(S)  /U(rHETA) 

Here  it  is  noted  that  there  are  no  pole  zero  pairs 
that  conveniently  cancel  so  that  the  transfer  function 
becomes: 

T3(S)  =  (1/IYY)  *DZIH*S*  (S  +  72.33)/ 

(S  +  35.1  +  j  27.84)*(s+  0.5085)*(S  +  0.36  +  J4.775) 

Making  the  simplification  used  in  th=  previous 
section,    the   result    is; 

T3  (5)     =     (72.33/IYY)  *DZra*S/(3**2    +    71.8*S    ♦    2055)* 
(S    +•    0.508)*(S**2    +•    0.72*S    +    22.74) 

Which  can  be  further  simplified  by  lookiag  at  the 
long  term  response  vice  initial  few  milli-seconds  to 
(S/(5    +   0.508))     =    2    so    that: 
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T3(S)  =  (142/IYY) *DZTH/(s**2  +  71.8*5  +  205S)*(5**2  + 
0.72*s  +  22. 74) 


4.   Tl  (S) 

Since  the  transfer  fandtions  are  the  same  except 
for  the  gains  the  solutions  are  the  same  except  for  the 
gains : 

T4<S)  =  (  1/HASS)  *DTHZ  *  S*(S  +  72.33)/DELTA 

where  the  same  simplifying  assumptions  may  be  applied. 

F.   SIMPLIFIED  BODE  PLOT  ANALYSIS 


*•   Simplifying  Assumption 
(see  fig.  12) 

The  assumption  is  made  that  the  craft  linearized 
frequency  response  transfer  functions  can  be  separated 
into  two  transfer  functions. 

The  first  block  G1  ( W, U, . .  . )  expresses  the 
disturbing  moments  ani  accelerations  generated  by  the 
craft  interacting  with  the  waves.  This  function  was 
developed  by  Thaler  and  Gerba  P  Uj  for  the  heave-only 
model.  In  that  simpler  model  the  function  31  was  a 
complex  function  of  encounter  frequency,  ingle  of 
incidence,  craft  length  compared  to  encounter  wavelength 
etc..  This  portion  of  the  overall  frequency  response  of 
this   model   is   left   for   future   work   due   to    time 
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constraints. 

The  second  block  32  (We)  represents  the  heave-pitch 
dynamics  of  the  model  as  developed  in  chapters  II,  III, 
and  IV,  which  respond  to  the  disturbing  moments  and 
accelerations  from  block  S1. 

The  Transfer  functions  T1   and   T2   (unsinplif ied) 
were  simulated  using  the  IBM-360  simulation  language  DSL 
with  S  =  jW  with  the  following  results. 

-•      S2^S  £i2i§  for  Pitch  Dynamics 

(see  figs.   15   and   16  ) 

Figure  15  represents  the  normalized  (numerator 
gain  constants  set  to  1)  gain  in  dB's  versus  Omega. 
Note  that  the  abscissa  scale  is  Omega  in  powers  of  10 
since  this  is  a  convenient  means  to  achieve  a 
logarithmic  scale  on  the  X-Y  plotter  output.  ^s  can  be 
seen  the  response  strongly  supports  the  simplification 
to  second-order  system.  These  plots  represent  the  full 
transfer  function  T1 (jW)  vice  the  simplified  system. 
Figure  16  further  supports  the  simplification  as  the 
phase  plot  starts  at  zeco  degrees  and  shows  only  a  very 
slight  positive  increase  prior  to  the  rapid  change 
toward  -180  characteristic  of  an  underdamped  system. 
Also  note  that  above  the  damped  natural  freguency  the 
nagnitude  plot  falls  at  -12  dB  per  octave. 

3.   Bode  Plots  for  Heave  Dynamics 

(see  figs .  17  and  1  3) 
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Figure  17  represents  the  normalized  magnitude 
versus  Omega  far  the  h  =  ave  dynamics.  The  system  is 
characteristic  of  an  overdamped  third  order  system  with 
small  gain  increases  at  the  craft's  natural  pitch 
frequency,  which  represents  the  pitch  dynamics  feeding 
back  into  heave  dynamics.  However  the  gain  magnitude 
curve  has  already  fallen  20  da  below  it's  lower 
frequency  response  gain  before  the  pitch  coupling  begins 
to  manifest  itself.  As  can  be  seen  in  fig.  18  the  phase 
plot  shows  the  same  third  order  pole  response  with 
slight  second  order  zero  influence  at  the  natural  pitch 
frequency. 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 


As  shown  in  Chapters  IV  and  V  the  simplified 
Pitch-Heave  model  developed  herein  is  a  very  reasonable 
simulation  of  the  XR-3  test  craft  for  the  given 
operating  conditions.  Heave,  Pressure,  and  C.G. 
Acceleration  were  in  vecy  close  agreement.  Pitoh  on  the 
other  hand  showed  some  difference  in  natural  frequency 
and  damping  factor.  It  is  felt  that  this  difference  is 
primarily  attributable  to  the  inclusion  of  a  redundant 
planing  force  calculation  in  the  6  DOF  model. 
Historically  the  results  of  the  6  DDF  computations  11 
for  pitch  response  have  been  very  similar  to  the  results 
for  the  linear  model,  in  both  natural  frequency  and 
damping  characteristics.  This  can  be  attributed  to  the 
fact  that  the  additional  planing  force  did  not  exist  in 
the  original  5  DOF  program. 

Since  one  of  the  major  concerns  with  the  6  DOF  model 
is  the  excessive  computational  times  involved  in 
sea-state  simulations  the  conclusion  may  be  drawn  that 
eventual  substitution  of  the  sidewall  and  seal  force 
linear  model  equation  in  the  5  DOF  program  will 
substantially  reduce  CPU  time,  with  some  resultant  loss 
in  computational  accuracy.  By  way  of  comparison, 
disregarding  compiling  time  and  output  time,  the  linear 
model  required  approximately  9  seconds  of  CPU  time  for 
the  5  second  simulation  while  the  6  DOF  model  required  9 
minutes  of  CP'J  for  the  5  second  run.  This  is  a 
computation  time  improvement  on  the  order  of  60  times. 

From  the  time-domain  plenum  gauge  pressure  data  it  can 
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be  seen  that  the  dynamic  variation  in  effective  plenum 
roof  area  has  a  noticeable  effect  on  the  plenum  pressure 
transient  behavior  as  compared  to  previous  results  14  . 
Ihe  initial  pressure  down  transient  was  faster,  dropped 
to  a  lower  value  ,  and  began  it's  recovery  more  quickly. 
Therefore  the  plenum  roof  dynamic  variations  strongly 
affects  the  C.G.  acceleration  characteristics  of  the 
craft  and  therefore  it's  transient  behavior. 

Ihe  dynamic  center  of  pressure  variation  developed 
herein  proved  to  be  a  very  draft  sensitive  phenomenon. 
Although  the  length  o£  the  other  lever  arms  as  modeled 
herein  change  by  a  similar  amount,  percentage  of 
variation  is  radically  different.  For  the  seals  the 
variation  is  on  the  order  of  inches  in  10  feet,  where 
the  variation  in  lever  arm  for  plenum  pressure  lift  is 
on  the  order  of  inches  in  a  single  foot  (or  less.)  This 
implies  that  dynamic  Z.  P.  variation  is  the  lominant 
factor  in  steady  state  pitch  angle  versus  steady  state 
draft. 

Planing  forces  as  included  in  this  model  are  now  in  a 
form  where  the  basic  factors  as  far  as  craft 
construction  characteristics  are  readily  recognized. 
This  force  provides  part  of  the  answer  to  craft 
pitch-down  characteristics  as  initially  investigated  by 
Reidel  C93  ,  with  the  remainder  of  the  moment  being 
developed  by  the  C.P.  variation  for  pressure.  The 
planing  force  provides  a  clue  as  to  how,  during 
construction,  the  craft's  pitch  amplitude  excursion 
characteristics  may  be  controlled. 

Several  areas  for  future  investigation  have  become 
apparent  in  this  woric. 

1.  Verification  of  the  effect   of   dynamic   center   of 
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pressure  variation  and  planing  forces  as  incorporated  in 
the  6  DOF  model  over  the  entire  speed  operating  range. 

2.  Incorporate  the  simple  modeling  for  sidewalls  and 
seals  into  the  6  DOF  model  and  investigate  the  trade-off 
between  CPU  time  and  conputat ional  accuracy. 

3.  A  closer  look  at  the  seal  forces  used  herein.  It 
is  possible  that  a  simple  gain  factor  needs  to  be 
included  in  these  equations  in  order  to  optimize  the 
steady-state  seal  forces  with  respect  to  actual  load 
data  as  measured  during  experimental  test  runs  of  the 
craft   10  . 

4.  Incorporate  sea  state  forces  and  moments  into  the 
simplified  model  by  simulation  of  the  forces  and  moments 
that  are  generated  by  craft  interaction  with  the  sea. 
Evaluate  the  accuracy  of  the  method  used  by  comparison 
with  data  generated  by  the  simple  model  and  the  6  DOF 
model  in  sea-state  conditions. 

5.  Analytical  investigation  of  the  craft  physical 
characteristics  that  have  predominant  influence  on  the 
natural  pitch  frequency  and  damping  by  using  the  gain 
terms  from  the  results  of  the  signal  flow  graph  analysis 
done  in  Chapter  V. 

6.  Use  the  simplifiai  model  to  develop  a  ride-comfort 
control  system.  lest  the  resulting  design  on  the  6  DOF 
model,  and  compare  the  results  obtained  with  those 
suggested  by  Rd'dR  C7]  .  It  has  been  known  ampirically 
for  some  time  that  the  craft  exhiaited  simple  second 
order  behavior  in  the  pitch  mode.  Mow  that  this  fact 
has  been  analytically  verified  the  implications  for  an 
automatic  control  system  for  pitch  response  are 
powerful.    A   simple   control   system  can  be  visualized 
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that  does  not  involve  either  complex  optimal  control 
system  design  work  or  hardware.  A  simple  acceleration 
and  rate  feedback  control  system  is  implied,  with  the 
possibility  of  adding  a  speed  input  to  the  system  to 
compensate  for  the  localized  linearity  of  the  craft 
characteristics. 
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APPENDIX  A 


LINEAR  SES  MODEL  USSR'S  GUIDE 


A.   General  Comments 

The  program  as  listed  following  Appendix  B  is  intended 
for  useage  via  the  hot  card  reader  at  NPS. 

It  is  intended  to  run  in  Fortran  G  with  Plotting 
option  (PROC  =  PORICLGP.) 

The  program  is  structured  on  the  assumption  that  the 
initial  and  final  state  variables  are  in  eguilibrium 
prior  to  data  entry,  but  will  run  regardless.  rha  force 
and  moment  residual  listings  will  indicate  the  value  of 
imbalance,  if  any. 

The   craft   physical   dimensions   are    alterable    by 
changing  the  cards  in  tne  craft  dimension  section  of  the 
main  program,  which  are   keyed   to   the   craft   geometry 
herein. 
B.   Data  Deck  Setup. 

The  /*  and  //GO.SYSIN  DD  *  cards  are  inserted 
following  the  and  card  ia  Subroutine  VERSAP,  foLlowed  by 
the  data  deck  as  described 

LPlot  option  card  (l?LOT)  ,  format  11.    IPLOr  =   1   = 
VESSATEC  X-Y  plotter  output. 
IPLOT  =  0  =  ?RT  PLOT  output 
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2.  List  option  card  (ILIST),  Format  11 

ILIST  =  1  =  Forces  and  moments, with  residuals  and  all 
sensitivity  coefficients  listed. 

ILIST  =  0  =  input  craft  operating  conditions  only 
listed . 

3.  Time  and  Integration  variables. 

TI,  TF,  NPTS,  NSTEP  Format  (2F1 0. H, 2110) 

Ti  =  problem  start  time  (normally  0.0) 
TF  =  Problem  stop  time  (seconds) 
NPTS  =  number  of  tabular  output  points  desired. 
NSTSP  =  number  of  integration  steps. 

note:  it  is  best  to  make  NSTEP  an  integer  multiple 
of  NPTS,  and  never  less  than  NPTS. 

4. Planing  coefficient  and  arbitrary  center  of  pressure 
input,  i.e.  HEAD  (  )  AKP,  PLCOEFF 
Format  =  2F10.4 

5.  Initial  Conditions  card. 

The  initial  craft  conditions  are  read  in  in  the 
following  order:  weight (Lbs .) ,  Draft  (in),  Pitch  angle 
(Deg.)),    PEBAR    (psf) ,    and    speed(fps.)       Format    is      5F10.4 

6.  Final  Conditions. 

The  final  conditions  ace  read  in  the  following  order 
and  format:  weight  (Lbs),  Draft  (in),  Pitch  angle 
(Deg),  and  PBBAR  (psf)  ,  Format  4F10.4.  Since  this  is  a 
constant  speed  model,  initial  and  final  speed  are 
assumed  agual. 

7.  Label  Cards. 

Label  cards  for  VERS&IEC  output  are  the  last  10  cards. 
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There  are  five  graphs  in  this  mode.  Draft(in)  ,  C.  G. 
Acceleration  (G*  s)  ,  Pitch  Angle  (Deg. )  ,  Pitch 
Acceleration  (Deg/SEC**2) ,  and  PBBAfi  (Each  graph 
reguires  two  label  cards,  Format  (6A8),  which  print  on 
two  lines  below  the  graph.  The  graphs  are  output  in  the 
order  listed  above. 
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APPENDIX  B 


DRAWINGS  AND  FIGURES 
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Figure  2    -       ASSUMED    GEOMETRY    <CONT> 
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Figure    4  -   PLENUM  CHAMBER  MO  DEL  IMG 
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figure    5  -   LINEAR  PLENUM  PRESSURE  TRANSIENT 
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Pigura   13  -   SIGNAL  FLOW  GRAPH 
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Figure   14  -   HEAVE  PITCH  TRANSFER  FONCTIDN 
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